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Abstract 



^ . A first-principles statistical theory is constructed for the evolution of two 

dimensional interfaces in Laplacian fields. The aim is to predict the pattern that 
J> ■ the growth evolves into, whether it becomes fractal and if so the characteristics 
of the fractal pattern. Using a time dependent map the growing region is confor- 
mally mapped onto the unit disk and the problem is converted to the dynamics 
O I of a many-body system. The evolution is argued to be Hamiltonian, and the 
^ I Hamiltonian is shown to be the conjugate function of the real potential field. 



Without surface effects the problem is ill-posed, but the Hamiltonian structure 
^ ! of the dynamics allows introduction of surface effects as a repulsive potential be- 



tween the particles and the interface. This further leads to a field representation 



^ , of the problem, where the field's vacuum harbours the zeros and the poles of 
CJ ■ the conformal map as particles and antiparticles. These can be excited from the 
vacuum either by fluctuations or by surface effects. Creation and annihilation of 
^ ' particles is shown to be consistent with the formalism and lead to tip-splitting 
■ and side-branching. The Hamiltonian further allows to make use of statistical 
mechanical tools to analyse the statistics of the many-body system. I outline the 
way to convert the distribution of the particles into the morphology of the inter- 
face. In particular, I relate the particles statistics to both the distributions of the 
curvature and the growth probability along the physical interface. If the pattern 
turns fractal the latter distribution gives rise to a multifractal spectrum, which 
can be explicitly calculated for a given particles distribution. A 'dilute boundary 
layer approximation' is discussed, which allows explicit calculations and shows 
emergence of an algebraically long tail in the curvature distribution which points 
to the onset of fractality in the evoloving pattern. 
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1. Introduction 

Notwithstanding the abundance of phenomenological knowledge the morpholo- 
gies of interfaces that grow in Laplacian fields are poorly understood theoretically. In 
many cases these interfaces exhibit a rich variety of convoluted patterns. Known ex- 
amples arc diffusion-limited aggregation (DLA), solidification of supercooled liquid, 
electro deposition and growth of bacterial colonies, to mention but a few. At present 
there is no sufficient theoretical understanding to allow for reliable predictions of 
the asymptotic patterns that such growths evolve into, starting from the basic equa- 
tions of motion (EOM). Existing analyses are either of an effective medium type [1] or 
employ renormalisation group techniques assuming similarity solutions and limit dis- 
tributions. Recently two of the latter generic approache managed to yield quite 
accurate values for the scaling of the growth size (or the mass) with the linear size 
of the growth[2][3]. 

A different direction to treat this problem was suggested exactly half a century 
ago[4]. It was proposed that two-dimensional free interfaces (i.e., in the absence of 
surface tension) be analysed by conformally mapping the growing region onto the 
unit disk and studying the evolution of the map rather than that of the interface. 
This idea was followed by Richardson [5] who discovered that such a map enjoys a 
set of independent constants of motion. Shraiman and Bensimon[6]took this issue a 
step further showing explicitly how the problem can be converted to a many-body 
system and pointing out that the free interface evolution is mathematically ill-posed. 
They demonstrated that the formalism breaks down after a finite time due to in- 
stabilities with respect to growth of perturbations along the interace on ever shorter 
lengthscales. Without the curbing effect of surface tension, irregularities develop into 
cusp singularities along the physical interface. There have been theoretical efforts to 
counteract this catastrophic sharpening by using small surface tension to cut off the 
short lengthscales in a renormalisable manner[7]. These, however, met with another 
difficulty: Such an ad-hoc inclusion turns out to constitute a singular perturbation to 
the EOM of the system. This, in turn, means that solutions based on this approach 
are strongly sensitive to initial conditions in the sense that very similar initial condi- 
tions can result in completely different morphologies [8]. Yet, numerous observations 
show that the final morpholgical properties are rather robust to the details of the 
growth and to changes in initial conditions. A different approach has been proposed 
recently to prevent formation of cusps by arguing that tip-splitting reduces high sur- 
face curvatures. Implementing this idea into the EOM of the equivalent many-body 
system results in production of particles [9], as will be elaborated on in this paper. 

A significant question in this context, that has not been addressed much, is 
whether the system supports a Hamiltonian or a Lyapunov function[10] [11]. This 
point is essential to a fundamental understanding of this highly nonequilibrium 
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growth process and therefore to the thrust of this paper, as will become clear be- 
low. Although, as formulated, the problem is known to support a set of conserved 
quantities[5][9][12], it is unclear whether these can assist in finding a Hamiltonian for 
the system. 

In this paper I first introduce the EOM of the interface. I show that the interface 
follows Hamiltonian dynamics and that the Hamiltonian is directly related to the 
physical two dimensional potential. I then derive the EOM of the equivalent many- 
body system, whose particles are the zeros and poles of the map. Using the existence 
of a Hamiltonian formulation, surface effects can be described as a field that either 
repels particles from the surface or gives rise to particles production. Both approaches 
overcome the ill-posedness of the problem and extend its range of validity to infinite 
time. The Hamiltonian also allows to analyse the statistical mechanics of the many- 
body system and in particular the spatial distribution of the particles inside the 
unit disk. I show how to extract information on the morphology of the interface 
from the distribution of poles and zeros. This includes a calculation of the moments 
of the growth probability distribution along the interface and the distribution of 
the curvature. I demonstrate the calculation for the case when the particles form a 
dilute gas near the unit circle and show that the distribution of the curvature along the 
interface develops an algebraic tail, implying non-negligible probabilities of formation 
of particularly sharp protrusions. This tail points to the onset of fractality and self- 
similarity, a feature that is not asumed a-priori. The third moment of the growth 
probability along the interface is calculated explicitly in terms of the distribution of 
the particles. This moment gives the fractal dimension of the pattern. 



2. The basic problem and the interface's EOM 

The fundamental problem of Laplacian growth can be formulated as follows. 
Consider a Jordan curve, 7(5), embedded in two dimensions which represents the 
physical interface. This curve is parametrised by < s < 27r, and is fixed at a given 
value of the potential field (electrostatic potential for electrodeposition or concentra- 
tion for diffusion controlled processes). A higher potential is assigned to a circular 
boundary whose radius, R, is very large compared to the growth size. The potential 
field, outside the area that is enclosed by 7 satisfies Laplace's Eq. 

V^$ = 0. (2.1) 

The interface is assumed to grow at a rate that is proportional to the local gradient 
of the field which is normal to the interface[13] 

= -V$ • h . (2.2) 
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This rate is assumed to be sufficiently slow so that at any time the Laplacian field can 
be considered to be static. Denoting hy ( = x + iy the physical (complex) plane, we 
now conformally map at each instant of time, t, the curve onto the unit circle (UC) 
in a mathematical z plane via ( = F{z, t). The time-dependent interface is recovered 
from the map by 7(s,t) = lim^^e" F{z,t). The field gradient along the curve is 
— V$(C) = —[d^{C)/dC]* = —i/{zF')*, where * stands for complex conjugate and 
the prime indicates derivative with respect to z. Using the fact that z — e'^^ = 1/z* 
on the interface, Shraiman and Bensimon [6] derived the EOM for 7: 

da{s,t) = -ids^is,t)[\ds^{s,t)\-'' + igis)^ . (2.3) 

The first term within the square brackets on the r.h.s. of (2.3) represents the normal 
growth rate as constituted by relation (2.2). The second term, however, is added by 
hand (but it is uniquely determined) and represents a tangential velocity, or 'sliding', 
of a point along the interface. Denoting the entire square brackets on the r.h.s. of 
Eq. (2.3) as the limit of an analytic function G{z, t) (whose explicit form is unique 
and explicitly determinable - see below), the EOM for F{z, t) becomes 

F = zF'G . (2.4) 

For reasons to become clear below let us also write down the logarithimc derivative 
of this equation with respect to z: 

The map needs in general to satisfy several constraints [9]: First, we want the topol- 
ogy of the boundary far away to remain unchanged under the map, which enforces 
lim^^oo F z. Second, the map must have no branch cuts. These two conditions 
are satisfied by the following general form 

where the time dependence appears in the scaling factor, A, and the locations of 
the zeros and poles of F'. The first constraint imposes 'charge neytrality', i.e., the 
number of poles should equal the number of zeros. The second constraint imposes a 
'dipolar neutrality', namely, '^^Zn = ^^Pn- It is possible to show that this form, 
and with a proper choice of the number N, enables to describe any initial simply 
connected cuspless curve that is allowed by this process. Therefore, this form is quite 
general and not merely a small class of maps[14]. Integrating (2.6) we have 



C = F{z,t) = A{t) 



N 



1 + ^RnHz-Pn) 



n=l 



(2.7) 
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The quantities Rn are the residues of the product in (2.6) at P„ (without the prefactor 
A), and 



n=l 



^ him -TT -Tm 

m=l ^"^"^ m=l 

By contour-integrating around the location of the poles and zeros in (2.4) and (2.5) 
one obtains their EOM: 

-A\t)z^ = zJgo + Y: + qJi - E ^ ^n'\{z}; {P}) 

m m 

-A\t)P^ = pJGo + Y: - fPiiZy, {P}) ■ 

(2.8) 

The behaviour of the particles under these these equations was discussed in detail 
by Blumenfeld and Ball[9]and will not be repeated here. From (2.4), and using the 
requirement that there are no branch cuts, one finds that the quantities A{t)Rn{t) are 
constants of the motion that are independent of each other and which are determined 
only by initial conditions. These constants are directly related to those found by 
Richardson [5] and mineev[12]. The time evolution of A{t) is straightforwardly found 
from (2.4) 

A{t) - A{t)Go . (2.9) 

This relation shows how this prefactor dependends on the distribution of the particles 
inside the unit disk. This rescaling factor is related directly to the fractal dimension 
of the pattern, if it becomes fractal, as will be discussed in section 5. 



3. The Hamiltonian structure 

Whenever a set of dynamical equations appears on the scene the first pertinent 
question is whether it supports a Lyapunov function or a Hamiltonian. If it does 
this gives access to a powerful bag of tools. I argue that the present probcess is 
indeed Hamiltonian and further relate it to the actual field The first hint that 
the process is Hamiltonian comes from the physical EOM (2.3), when only normal 
growth is considered: 

S s 

dais,t) = -idsi{s,t) \dsi{s,t)\~'^ = -^T— • (3.1) 
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This expression and its complex conjugate are equivalent to Hamilton's equations. 
This relation suggests that the (rescaled) length of the actual growth, s, may play 
a role of a Lyapunov functional, while 7 plays the role of a field, whose real and 
imaginary parts are the canonical varialbles. Although Eq. (3.1) is not formally 
correct this appealing interpretation begs the question whether the EOM for the 
map, F, also supports such a structure. It is not difficult to show[15]that this is 
indeed so: Multiply the complex conjugate of (2.3) by and rewrite its imaginary 
part in terms of the map [4] 

Im{F;F,} = l . 

Writing F = 00 + ix, with u and x real functions, this equation asserts that the 
Jacobian of the transformation from the coordinates t — s to the coordinates — x 
is unity, d {uj,x) /d {t, s) = 1. It is then straightforward to show that when z ^ 



ds . ds ^ ■ . ds .„ 

Since s = Im{ln2;} Eq. (3.2) generalises to 

where $ is the complex potential in the physical plane. Thus the trajectories of the 
system follow the stream lines, which are the conjugate of the equipotential lines in 
the physical plane. 

Yet another way to obtain a similar structure is as follows: Define the complex 
function ^' = F' + iz. By manipulating Eq. (2.5)[16]one obtains 

.dHo 

-dt = -'d^ ' ^^-"^ 

where Hq = zF'G, which appears on the r.h.s. of (2.5) is complex and therefore does 
not enjoy a convenient translation into Hamilton's equations. 

The main point regarding the above arguments is not as much the exact form 
of the Hamiltonian but rather that the p.d.e. that governs the interface's evolu- 
tion indeed follows Hamiltonian dynamics. In other words, given an initial value 
of the Hamiltonian (the 'energy') the system then follows a trajectory that keeps 
this value constant. Since the many-body formulation is an equivalent description of 
the growth process it follows that the system of poles and zeros must also keep this 
quantity constant and hence the latter also supports a Hamiltonian structure. To 
find the many-body Hamiltonian one inverts relation (2.7) to express In 2; in terms 
of F and then Im{ln2} = ^{F, F*). I will not pursue this direction further here but 
rather argue that the very existence of a Hamiltonian already paves the way to much 
progress. 
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4. Surface effects 

Turning to consider surface effects, it has already been mentioned that without 
surface tension (or capilary forces) cusps form along the interface due to instability of 
small corrugations [17]. But real growth processes clearly do not admit cusps. In the 
cases that concern us here this is because the system has to expend a macroscopic 
surface energy as the curvature increases. What does the procedure of cusp formation 
correspond to in the many-body system? A local protrusion (the incipient cusp) is 
caused, in the many-body system, by a zero approaching the UC. Thus prohibition 
of high curvatures naturally corresponds to keeping that zero from approaching the 
UC too closeely and hence to an effective repulsive potential between the particles 
and the interface defined by the UC. It should be stressed that only the existence of 
a Hamiltonian makes it possible to use the term 'repulsive potential' with any proper 
meaning. 

There are different ways to incorporate this idea into the theory: One is by 
introducing a surface potential term in the many-body Hamiltonian[16]. Although 
this may sound somewhat difficult since we don't know the exact Hamiltonian of 
the many-body system, one can nevertheless insert such a term in the p.d.e. (3.4), 
H = Hq + V, and derive the modified EOM for the particles. The choice of the surface 
potential term determines the nature of the growth to a large extent. The stronger 
the repulsion, the smoother the resulting interfaces. An example of a possible simple 
repulsive potential term is 

V = a lim ln[K{{Z},{P})] , (4.1) 

where K is the (complex) curvature in terms of the locations of the zeros and the 
poles [9] 

^7 P . 

K(., izy, {P}) = ^H^^ + El-^ -T^J}' (^-2) 

whose real part is the physical curvature. This particular surface potential is simple 
in that it contributes a constant repelling term in the EOM of the particles, (2.8). It 
has recently been shown that a term that diverges as a particle approaches the UC 
would do better to describe the physics[16]. 

Another approach is to make a deeper use of the fact that the particles move 
in fact in a field. The field can have a vacuum that can accommodate particles 
and antiparticles (zeros and poles). 'Exciting' the vacuum (say, by fiuctuations) can 
then effect creation and annihilation of zeros and poles, a mechanism that allows 
for tip-splitting and side-branching [9]. But before considering such a far fetching 
interpretation we need to convince ourselves that such a picture is consistent with 
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the present formalism and that it does not contradict any of the basic premises. To 
this end consider again the derivative of the map (Eq. (2.6)) and rewrite it at the 
initial moment (say, to) in the, seemingly redundant, form 

The second product is unity at to and at any time thereafter under the free interface 
evolution because all its terms simply cancel out. It is straightforward to see that, 
under the free-growth EOM, at any time t > to F' will retain this form with only the 
position of the particles in the first product and the value of A changing. Consider 
now what happens when a zero, Zi, and a pole, Pj, collide. At the instant of collision 
the particles occupy the same location and therefore their corresponding terms, (z — 
Zi)/{z — Pj), cancel out in the first product on the r.h.s. of (4.3). This is exactly 
an annihilation event. Such an event conserves the balance between the numbers of 
zeros and poles, so that charge neutrality is not violated. It also does not violate the 
dipolar neutrality, as can be immediately verified. Turning to production events, the 
second product on the r.h.s. of (4.3) can now be interpreted as a 'vacuum' of pairs of 
zeros and poles that do not manifest unless they are 'excited'. One can envisage two 
routes for this to occur: i) A fiuctuation, of whatever origin, can virtually separate 
such a pair; ii) A particle that moves with a high kinetick energy can knock the pair 
apart. Particles with high velocities are usually zeros that are close to the unit circle 
and therefore give rise to locally high curvatures. 

So the form (2.6) is also naturally suited for creation. Immediately after a zero- 
pole pair has been excited their locations are very close and it is quite straightforward 
to verify from the EOM that such a close pair interact repulsively, pushing apart and 
hence maintaining their identities. Differently expressed, once created the particles 
are stable. Thus a pair, say at Vj, can be excited into two individual particles at 
Zj = Tj + Si and Pj = Tj + 62- Under production, as under annihilation, the 
basic constraints need to be satisfied. Namely, an excitation is in pairs for charge 
neutrality, and dipolar neutrality is staisfied by imposing a relation between the 
locations of the excited particles (I should mention that excitation by field fluctuations 
has to occur in quartets rather than in pairs due to the dipolar constraint). In a 
particular implementation of this idea Blumcnfeld and Ball[9]proposed that excitation 
of a new pair is triggered by the proximity of a zero, say Z^, to the unit circle 
(scenario ii above). Such proximity leads to a high curvature in front of Zn at 
s = a,rg{Zn} and the new zero-pole pair is excited once the local curvature reaches 
a threshold value. To satisfy dipolar neutrality the new pole occupies the location 
of Zn prior to the production event, while the other two zeros are equidistantly and 
oppositely situated around the pole. It turns out that the orientation of the zeros 
is not limited by the constraints on the system and needs to be imposed following 
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another criterion. Blumenfeld and Ball introduced an energetic criterion along the 
following lines: Recognising that surface energy increases with increasing curvature, 
the location of the particles needs to minimise the local curvatute at s. As it happens, 
this minimum corresponds to placing the two zeros in the azimuthal direction about 
the location of the newly born pole. In the physical plane, such an event constitutes 
exactly tip-splitting. Placing the zeros along a radial ray equidistantly from the pole 
maximises the local energy and corresponds to side-branching in the physical plane. 
Since the system would rather minimise its local energy, these results suggest that tip- 
splitting may have an energetic, rather than only stochastic, origin. This particular 
procedure is not universal and other systems may follow other criteria for particles 
production. It should be emphasised that it is quite plausible that fluctuations of 
the field superimose on this mechanism and stochastically induce tip-splitting and 
side-branching. Further studies in this direction are currently being carried out and 
will not be elaborated on here. 

5. Noise and statistical analysis 

The foregoing assumed mostly a deterministic evolution. It is known[18]that 
the dynamical EOM of the p.d.e. leads to a chaotic growth in the sense that if 
one starts from very close initial conditions, one ends up very quickly with different 
deterministic structure. This implies an efficient spread of the solutions in phase space 
which, combined with the good fortune of having a Hamiltonian dynamics, ensures 
the existence of a Gibbs measure. Namely, one expects to be able to construct a 
partition function 



where VTZ = fln^i (PZ^d^Pr, is a an infinitesimal volume in phase space and the 
Hamiltonian is that of the many-body system. The quantity (3 reflects the 'noise' in 
the system and is obtained as a Lagrange multiplier by imposing an average 'energy' 
constraint on the distribution. Expectation values of quantities such as the moments 
of the curvature and moments of the field gradient, |V$|, can now be found via 



Since, at present, the explicit form of H is unknown, this approach is not easy to 
implement. 

Alternatively, we can construct a master equation for the evolution of the dis- 
tribution of the particles, N{{Z}, {P})i using Liouville's theorem 




(5.1) 




(5.2) 




(5.3) 



9 



where F represents collisions and noise. To use Liouville's theorem let me confine 
myself here to systems with conserved number of singularities (no particles produc- 
tion). An extension to a nonconserved number of particles is not difficult once a 
self-consistent renormalisation is introduced and will be discussed elsewhere. Ex- 
pecting a steady state distribution after rescaling the growth by A{t) we can discard 
the explicit time derivative. A solution of this equation yields everything there is 
to know about the interface's statistics. Unfortunately, as common in statistical 
mechanics, a general exact solution is impossible. Rather than trying to solve this 
equation in some approximation, let me demonstrate how such a solution can be 
converted into information about the physical interface. 

First, the distribution of the curvature, Vk, can be derived from J\f using the 
relation 



Vk^ /Ar({Z},{P}) Ji^i-^ 



TV 



1 + Re 



1 1 



vn 



(5.4) 

where Ki is the measurable curvature along the interface and 5 denotes Dirac's 
delta-function. More generally, the distribution of any morphology-related quantity, 
X, V{X), that is expressible in terms of the locations of the particles can be found in 
this manner. As an explicit example of using the statistics, I now turn to calculate 
the values of the moments of the growth probability distribution along the interface. 
This calculation has been recently carried out by Blumenfeld and Ball[9]and is only 
briefly reviewed here. These moments play a central role in pattern formation and 
growth, mostly because they were shown to lead to an asymptotically stable mul- 
tifractal function that is independent of initial conditions or details of the growth. 
The probability that growth occurs at a point s along the interface at some time t is 
proportional to the local gradient of the field 

p{s) = Co\W^s,t)\= lira Co \F'{z,t)\ ; Cq = 1/1 \V^{l)\dl . 



Therefore the quenched moments of this distribution are 



^ (b = — (f |F'(^)|'"^— . (5.5) 



'7 ^'^''J\z\=l 

Substituting from Eq. (2.6) we have 



M, 



n=l ^ n=l 



The second product in the integrand, = J{z)^^ i)/'^^ is analytic within the unit disk 
while the first contains N poles of order {1 — q)/2. For q > 1 these poles are located 
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at the zeros of the map, while for < 1 the poles of the integrand coincide with the 
poles of the map. In both regimes a simple pole also exists at the origin. Thus it is 
straightforward to evaluate this integral for odd values of q: 



M. 



1 



N 



1 



TV 



A{t) 



2u 



zJ{zY 



n 



Pk 



(5.7) 

where i^=(Q — l)/2isa positive integer number. Other values of q > 1 can be ob- 
tained either by taking into account explicitly the branch-cuts that are involved in the 
calculations, or by interpolating between the integer moments. However, it is known 
that the knowledge of all the odd moments of the probability density of a measure on 
a finite support is sufficient to determine that probability density uniquely [19], and 
therefore the present calculation suffices. Of particular interest is the third moment 

N ^ N 



Ms = A{ty 



Go/A\t) , 



(5.8) 



n=l 

where the last step makes use of the definition of Gq in section 2. Using Eq. (2.9) 
we then obtain the exact relation between Ms and A{t), 

Ms = A{t)/A%t) . 

The third moment is directly related to the fractal dimension[20][21]by 

Df = In Ms/ In Rg , 

where Rg is the linear size (the radius of gyration or the radius of the equivalent 
circular capacitor) of the growth. Df can be evaluated once the distribution of the 
particles is known. Moreover, if indeed the pattern turns fractal then this quantity 
should asymptope to a pure number. This is another manifestation of the first- 
principles nature of the present approach: It is the result that tells us whether the 
morphology becomes fractal without having to assume such a solution from the out- 
set. 

Calculating for odd values of q < 1 leads to an expression similar to (5.7). Since 
all these quantities depend on the locations of the particles one can find their distri- 
bution over many growth realisations via an integral similar to (5.4). An interesting 
observation is the following: Comparing the expressions for the positive and negative 
moments shows that these enjoy exactly the same form with the location of the zeros 
interchanged with the locations of the poles. This, combined with measurements on 
DLA that show a distinct difference between the negative and positive moments of 
the growth probability distribution implies qualitatively different spatial distribution 
of the two species of particles. A partial confirmation of this conclusion can be indeed 
observed in numerical calculations where the trajectories of the poles and the zeros 
display markedly different behaviours [9]. 
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6. The dilute boundary layer approximation 

To gain insight into some features of the morphology let me now discuss briefly 
the distribution of the curvatures using a dilute boundary layer approximation for 
the distribution of the particles. The motivation behind this approximation is as 
follows: Prom Eq. (4.2) one can observe that the curvature is dominated by the zeros 
that are closest to the UC. Namely, considering a zero, Z„ = (1 — p)e'*" with p -C 1, 
the curvature in front of this zero is (see Eq. (4.2)) 

i^l(Sn)~C/p' , (6.1) 

where, to a good approximation, C is independent of cither p or the locations of the 
other particles. A similar consideration for a pole near the UC shows that the local 
curvature is regular in p. We can therefore neglect the contribution of poles in the 
following approximation. Since these are the zeros with p <^ 1 that dominate the 
curvature we consider only the zeros within a ring 1 — pc < \A < 1 and assume that 
their density is dilute and isotropic (the isotropy assumption can in fact pertain only 
to a discrete number of global arms). It can be easily shown that this corresponds 
to analysing the exposed parts of the physical growth. By assumption then the 
curvature at Sn depends mainly on the radial location of whose distribution is 
Vnip). Thus 

Vk = Vn{p)dp/dK = Const. Kf'Vn ((C/Ki)^/^) . (6.2) 

In the absence of a solution to the master equation (5.3) we have no information 
on Vn- Nevertheless, we can observe that even if this distribution is well behaved 
(i.e., all its relevant moments are finite) the distribution of Ki exhibits algebraically 
decreasing tails for high curvatures. This immediately suggests the onset of fractality, 
which generically originate from such tails. It is important to note that nowhere 
along the construction of the theory did we assume fractality or self-similarity. Yet, 
an algebraic tail appears which can easily generate such a structure. The moments 
of K are dominated by the long tail and satisfy Hq ~ Kmax"^ ~ a^/^~^, where a is a 
small cutoff lengthscale. With growth the cutoff radius reduces by a factor of XjAi^) 
and therefore p^ ~ A(t)'^~^/^. It is quite plausible that Vn also introduces algebraic 
tails, directly affecting the behaviour ofVk- A more accurate analysis is only possible 
once we have a solution to Eq. (5.3) and attempts in this direction are being carried 
out. 

Before concluding this section I should mention that it is possible, using this 
approximation, to also evaluate the aforementioned moments of the growth proba- 
bility, Mg[9]. Such an evaluation shows that these moments depend strongly on the 
negative moments of the distribution of zeros' distances from the UC. If indeed the 
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patterns becomes multifractal, such an evaluation should yield the left hand side of 
the multifractal spectrum (the regime governed by the exposed parts of the growth) 
directly. 

7. Discussion and concluding remarks 

To conclude, a first-principles theory for two dimensional Laplacian growth was 
described. The growth was shown to be Hamiltonian with trajectories of constant 
energy corresponding to the stream lines in the two-dimensional physical field. This 
was argued to indicate that the equivalent many-body system is also Hamiltonian. 
Although not discussed here, an extension of this formalism to a continuous density 
of zeros and poles is possible and has been carried out.[16]The resulting nonlocal dy- 
namical equations, however, currently seem too complex for an analytical treatment. 
Surface energy was shown to lead to a field that acts on the particles in the mathe- 
matical plane. The field concept could be incorporated either as a repulsive potential 
between the particles and the unit circle or by giving rise to particles production. 
The latter corresponds to tip-splitting and side-branching. Thus this seems a nat- 
ural mechanism to prevent cusp formation along the physical interface and connect 
surface efi^ects to tip-splitting. I discussed the distribution of the zeros and poles 
inside the unit disk and related it explicitly to the morphology of the interface. The 
main advantage of this approach is that the nonequilibrium growth is describable in 
statistical mechanics formalism, the language of equilibrium phenomena. The dis- 
tribution of the curvature was analysed and an exact calculation of the moments of 
the growth probability distribution along the interface was discussed. To address the 
fundamental issue of onset of scale-invariance and self-similarity, I used an approxi- 
mation which showed that the curvature distribution develops an algebraic tail, which 
naturally gives rise to a fractal pattern. It is this author's belief that this approach 
and its possible generalisations are very promising as a framework for constructing 
theories for the morphology of growth processes in other dimensionalities and in fields 
that satisfy equations other than Laplace's. 
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